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Abstract. We consider a spatially homogeneous and isotropic system of Dirac 
particles coupled to classical gravity. The dust and radiation dominated closed 
Friedmann-Robertson- Walker space-times are recovered as limiting cases. We find a 
mechanism where quantum oscillations of the Dirac wave functions can prevent the 
formation of the big bang or big crunch singularity. Thus before the big crunch, the 
collapse of the universe is stopped by quantum effects and reversed to an expansion, 
so that the universe opens up entering a new era of classical behavior. 

Numerical examples of such space-times are given, and the dependence on various 
parameters is discussed. Generically, one has a collapse after a finite number of 
cycles. By fine-tuning the parameters we construct an example of a space-time 
which satisfies the dominant energy condition and is time-periodic, thus running 
through an infinite number of contraction and expansion cycles. 



1. Introduction, the Einstein-Dirac System 

Near the big bang or big crunch singularity, the matter density and the space-time 
curvature become arbitrarily large. It is generally believed that in this regime, quan- 
tum effects should come into play, which could even prevent the formation of the 
singularity. This effect has first been analzyed in [7] for classical gravity coupled to 
a second quantized matter field; see also [9]. For a quantized gravitational field, this 
effect has been studied in [6], and it was worked out in more detail in the framework 
of string cosmology ^Oj and in loop quantum gravity [IJ. In contrast to the above 
approaches, we are here more modest and work with classical gravity coupled to Dirac 
wave functions, without using second-quantized fields. We find a new mechanism, 
based on the oscillations of the spin of the matter field, which tends to prevent the 
formation of a space-time singularity. The mechanism can be understood as a spin 
condensation effect. An advantage of our model is that it is very simple and can be 
analyzed without any approximations. The Dirac field could be interpreted physi- 
cally as a new fermionic particle which becomes significant on the cosmological scale. 
Alternatively, the Dirac field can be understood similar to a "quintessence" as 
effectively describing unknown fields acting in our universe on the large. Describing 
such a cosmological field by the Dirac equation seems physically natural and changes 
the behavior of the cosmological model drastically near the big bang or big crunch, 
where quantum oscillations can prevent the formation of space- time singularities de- 
pending on the initial conditions. Such a bounce back may happen several times, 
before, generically, the solution collapses. For brevity we here restrict attention to the 
case of a closed universe. However, the flat and hyperbolic cases are less interesting, 
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since they expand infinitely. We remark that the hyperbolic case allows for solutions 
with negative energy density where the matter turns to be repulsive (see [1]). 

For the derivation of the equations we follow the standard approach in [HI [3] where 
we work with a torsion free connection (for the effect of torsion see the cosmological 
model in Einstein-Cartan theory p^J). The Einstein-Dirac (ED) equations read 

R] - ^R 6} = Sttk T] , (P - m)^ = , (1) 

where is the energy-momentum tensor of the Dirac particles, k is the gravitational 
constant, P denotes the Dirac operator, and ^ is the Dirac wave function (we always 
work in natural units h = c = 1). For the metric we take the ansatz of the closed 
Priedmann-Robertson- Walker (FRW) geometry 

ds'^ = dt^ - R^{t)da'^, 

where t is the time for an observer at rest, R is the scale function, and da is the line 
element on the unit 3-sphere, 

da^ = + r^d^^ + r2 sin^ ^ dip' , 

1 — 

where r G (—1,1) is a radial variable, and {'&,ip) € (0, §) x [0,27r) are the angular 
variables. The Dirac operator in this metric can be written as (see [3]) 

+ 2R{t) ) + Rit) V-^S3 

where = ( J -i ) is the usual 4x4 Dirac matrix, and P53 is the Dirac operator 
on the unit 3-sphere. The operator Vgs has discrete eigenvalues A = ±|,ib|,..., 
corresponding to a quantization of the possible momenta of the Dirac particles. For a 
given eigenvector -0;^, we can separate the Dirac equation with the ansatz 



= R{t)-'^ 



SlTK 



thus describing the time dependence of the wave function by two complex functions a 
and f3. In order to be consistent with the homogeneous and isotropic ansatz of the 
metric, the Dirac spinors must also be in a homogeneous and isotropic configuration. 
In analogy to the method employed in [5] for spherical symmetry, this can be achieved 
by taking an anti-symmetrized product of the wave functions ([3]), where ipx runs over 
an orthonormal basis of the eigenspace corresponding to a fixed A, where the symmetry 
is expressed by the fact that J2 |V'a('"5 "f?, 0)P = A^ — with the sum running over all 
corresponding angular quantum numbers. We thus obtain a fermionic Hartree-Fock 
state composed of A^ — | particles. The energy momentum tensor Tj can be derived 
in a similar way as in [5], see Equ. (4.4). By spherical symmetry the components 
T*,r^,r^ have to vanish as well as the components T^,T^,T^, which is due to the 
homogeneity in space. Therefore Tj is a diagonal matrix where the space components 
are equal, again due to the fact that our system is homogeneous in space. Analogously 
to [Sj we derive the time component of the energy momentum tensor, which reads 



T,' = ET^c.{n»(a, + ||)l*.)(*.l} 
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where the sum runs over the angular momentum quantum numbers, whose indices are 
suppressed for notational reason. Using ([3]) and the fact that the sum of \ipx{r, is 
constant in space, since the sum runs over the angular quantum numbers corresponding 
to A, we obtain that 



2A — 
m[\a\'^ - \P\'^] - — Re(a/3) 



(4) 



Substituting the ansatz ([3]) and (jl]) into the ED equations ([T|) , we obtain the following 
system of ODEs in the spinors (a, /3) and the scale function R (for a detailed derivation 
see i4j), 



'dt UJ " l-A/i? -m ) U ) 



^' + 1 = + «/?)• (6) 



Let us mention that the equations obtained from the spacial components of the energy- 
momentum tensor are automatically satisfied due to the continuity equation, for which 
reason we omitted its exact form. Let us emphasize again that all particles in our model 
have the same momentum A and therefore give rise to the same spinor equation ([5]), 
meaning they are represented by the same spinor (a, f3). This is what we mean by the 
above-mentioned notion of spin condensation. 

Note that ([6]) determines R only up to a sign, and at first sight this seems to lead to 
an ambiguity whenever R becomes zero. However, R becomes uniquely determined by 
demanding that be twice differentiable, implying that R must change sign at every 
zero of R. Normalizing the probability integral to one, we obtain the condition 

\a? + = ^ (a^ - ^) . (7) 

Differentiating and substituting (l5|), one sees that this normalization condition is in- 
deed time independent. 

A short calculation shows that the equations (O [6|) are for any p > invariant under 
the scalings 

in 

R — > pR , t ^ pt , m — > — , A — > A , (a, /3) — > /> (a, /3) . 

P 

Since this scaling changes the norm of the spinors, it is no loss of generality to re- 
place d?]) by the simpler condition 

\af + \P\^ = l. (8) 

An other way of understanding the above rescaling is that we choose the units for the 
gravitational constant such that the right side of d?]) equals one. 

The resulting ED system involves the physical parameters m (the rest mass of the 
Dirac particles) and A (the Dirac momentum, also related to the number of particles), 
as well as two free parameters to set the initial conditions of the spinors a and /?. 
These four parameters describe all the physical configurations of our system. 
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2. The Bloch Representation 

For the analysis of the system of ODEs (O [6]) it is convenient to regard the spinor 
(a, P) as a two-level quantum state, and to represent it by a Bloch vector v. More 
precisely, introducing the 3- vectors 

(where a are the Pauli matrices, and 61,62,63 are the standard basis vectors in R^), 
the ED equations become 

ij = bAv, i?2 + i = —b-v 

2R 

(where 'A' and '•' denote the cross product and the scalar product in Euclidean R^, 
respectively). To further simplify the equations, we introduce a rotation U around the 
62-axis, such that b becomes parallel to ei, 

2 



Ub = - + m2i?2 , 
R 

Then the vector w := Uv satisfies the equations 

^ ' (9) 



w = dAw , if + 1 = — 3_ y^A^ + m2i?2 ^ 



where 



Ami? R 



d:=- VA2 + m?R^ 61 - ^ - 62 . 10 

R A2 + m?R^ R ^ ' 

We refer to the two equations in Q as the Dirac and Einstein equations in the Bloch 

representation, respectively. Notice that equation (|1U|) can only have solutions if wi is 

negative. 

In the Bloch representation, one can most easily recover the dust and radiation 
dominated FRW-geometries. In the two limiting cases m ^ and A ^ 0, the second 
term in (jlOp drops out, and thus the ED equations simplify respectively to 

w = J— J- 61 A tZ; , R^ + 1 = — 75I wi (radiation dominated universe) (11) 

R R 

w = 2m ei Aw , R + 1 = —— wi (dust universe) . (12) 

R 

In both cases, the Dirac equation describes a rotation around the fixed axis ei. This 
implies that wi is constant in time, and thus the Einstein equations reduce to the 
well-known Priedmann equations for a radiation dominated universe and a dust uni- 
verse, respectively. Note that the components W2 and do oscillate due to quantum 
effects. These oscillations can be understood as the familiar "Zitterbewegung" of Dirac 
particles. But these quantum oscillations do not enter the Einstein equations. 

Away from the above limiting cases, the function wi is in general not constant, 
but it takes part in the oscillations as described by the first equation in ([9]). As a 
consequence, the dynamics can no longer be described by a single Friedmann equation, 
but only by the coupled system of Einstein-Dirac equations ([9]) . The relevant length 
scale is characterized by the radius 

i?qu = - • (13) 

m 
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Figure 1. Quantum oscillations preventing the big crunch for A = |, 
m = 10.5448 and wi(i?max) = -0.2675. 



On this length scale, the quantum oscillations do enter the Einstein equations, leading 
to effects which go beyond the scope of classical cosmology. We point out that the 
radius -Rqu can be much larger than the Planck length. Thus, in contrast to [H llOj . 
our mechanism is not related to quantum gravity; instead we are working consistently 
within the framework of Einstein's general theory of relativity. 

3. Quantum Oscillations Preventing the Big Crunch 

Using a standard ODE solver, we shoot for numerical solutions starting at the point 
where the scale function R reaches its maximum i?max) solving forward and backwards 
in time, until we reach a big crunch or big bang singularity, respectively. In Figure [T] 
a typical solution is shown. The function R increases after the big bang similar as in 
the classical FRW solution up to its maximal value Rmax, where it starts decreasing. 
The classical contraction stops on the scale R ~ -Rqu; where the quantum oscillations 
become important. This leads to an oscillatory behavior of R. As a consequence, 
the sign of R may flip at some minimal radius i?min; where the universe again begins 
to expand, going over to a classical FRW-like space-time. The functions W2 and ws 
oscillate, with slightly varying amplitude. The function wi is nearly constant in the 
classical region, but becomes oscillatory near the quantum regime R ~ Rq^ (see Fig- 
ure [T]) ; its nonlinear interplay with the scale function R leads to the surprising effect 
that the big crunch is avoided. 

The behavior in the quantum regime depends crucially on the phase (j) of the Dirac 
oscillations, defined by = arctan(t(;2/ti'3). In Figure [2] the minimal radius -Rmin is 
plotted as a function of the initial value of (j). If i?min > 0, the big crunch is avoided 
and R "bounces back," whereas in the case i?min = the big crunch appears despite 
the quantum effects. Considering the phase (p as unknown, one can give the result of 
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Figure 2. Dependence of the minimal radius on the initial angle (p 
for A = |, m = 10.5448 and t(;i(i?max) = -0.2675. 
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Figure 3. Quantum oscillations preventing the big crunch for A = |, 
m = 10.5448 and wi{Rm^y,) = -0.0608. 

Figure [2] a statistical interpretation: For a random initial phase, this bounce appears 
with a finite probability (in the example of Figure [21 this probability has the value 
p = 0.276224). 

Let us briefly discuss the qualitative dependence of the bounce on the free parame- 
ters i?max! -Rqu) A, and (j) of our model (for a detailed analysis we refer to [4J). First of 
all, the probability of preventing the crunch can be increased by choosing |rf i(i?max)| 
smaller. This is shown in the example Figure [3l where the bounce occurs even with 
probability one. We point out that small |T«i(i?max)| does not imply that i?max is also 
small. Namely, as one sees from the Einstein equation in ([9]) by setting R equal to zero, 
a small value of |?i'i(i2max)| can be compensated by choosing m or A large. This leaves 
us the freedom to choose -Rqu, (fT3l) . at will. In this way, one can construct space-times 
with a high probability of bouncing in the physically relevant case where -Rmax is the 
size of our universe and i?qu is a microscopic length scale. 

An interesting limiting case of our ED equations is obtained by dropping the first 
summand in (jlOp . This is a good approximation in the so-called scale dominated 
region where R decreases so fast that the oscillations of W2 and w-^ can be neglected. In 
particular, the scale dominated region describes the bounce in the physically interesting 
limit i?max oo (for fixed A). In this limiting case, we can integrate the Dirac equation 
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Figure 4. A time-periodic ED cosmology for A = |, 
t^lj^o = (0.690234,0,-0.723586). 



m = 21.3286, 



in (0) explicitly and express w as a function of R, 

w{t) = 




where ^{t) = — arctan[i?(t)/i?qu] + arctan[i?(to)/-Rqu]) and to is the time from which 
on the scale dominated approximation applies. For the probability p that the crunch 
is avoided one obtains 

1 / kl(^max)| Rqu \ 

p = — arccos — , — — ^ — . 

^ Vv'l-W^l(i?max)2 ^max ) 

This probability tends to 1/2 as -Rmax tends to infinity, showing that the bounce is of 
significance even for fixed |'Wi(i?niax)| in an arbitrarily large universe. 

4. Time-Periodic Solutions 

By iteratively adjusting the starting values at t = 0, we constructed time-periodic 
solutions. In Figure H] an example of a time-periodic solution is shown. Except near 
the "quantum turning points" between collapse and expansion, the space-time is well- 
approximated by the classical dust-dominated FRW universe. It is remarkable that 
the dominant energy condition is fulfilled even near the quantum turning points, as 
is shown on the right of Figure [H We point out that the periodic solutions were 
constructed by fine-tuning the initial condition and are thus not generic. A generic 
solution collapses after a finite number of cycles. 



5. Conclusion 

A new class of solutions of the Einstein equations coupled to Dirac spinors was 
constructed. These solutions satisfy the dominant energy condition and are thus phys- 
ically relevant. Our model reveals in a simple and explicit setting a general mechanism 
which tends to avoid space-time singularities, such as the big bang or the big crunch, 
if the quantum mechanical nature of matter is taken into account. By fine-tuning the 
initial conditions our model even allows for periodic solutions, an eternal universe with 
an infinite number of cycles. 
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